LOCAL BROWNIAN PROPERTY OF THE NARROW WEDGE 
SOLUTION OF THE KPZ EQUATION 



JEREMY QUASTEL AND DANIEL REMENIK 



Abstract. Let 'H{t, x) be the Hopf-Cole solution at time t of the Kardar-Parisi-Zhang 
(KPZ) equation starting with narrow wedge initial condition, i.e. the logarithm of the 
solution of the multiplicative stochastic heat equation starting from a Dirac delta. 
Also let W^{t,x) be the solution at time t of the KPZ equation with the same noise, 
but with initial condition given by a standard two-sided Brownian motion, so that 
'hf"^{t, x)—'W^{Q, x) is itself distributed as a standard two-sided Brownian motion. We 
provide a simple proof of the following fact: for fixed t, 'H{t, x) — {T-C'^{t, x)—'hC^[t, 0)) 
is locally of finite variation. Using the same ideas we also show that if the KPZ 
equation is started with a two-sided Brownian motion plus a Lipschitz function then 
the solution stays in this class for all time. 



1. Introduction and statement of the results 
The KPZ equation 

(1.1) dt'H = -]^{d^,nf + ]^dln + w, 

was introduced by Kardar, Parisi, and Zhang [9] as a model of randomly growing inter- 
faces. Here 'W{t,x) is Gaussian space-time white noise, ¥.{W{t,x)W{s,y)^ = 6s=tSx=y 
(see Section 1.4 of [2] for a precise definition). It is expected that the one-dimensional 
KPZ equation appears as the weak asymptotic limit of a large class of stochastic inter- 
acting particle systems/growth models, including directed random polymers, stochastic 
Hamilton-Jacobi-Bellman equations, stochastically perturbed reaction-diff'usion equa- 
tions, stochastic Burgers equations and interacting particle models, and it is in fact 
rigourously known to describe the fluctuations in weakly asymmetric exclusion pro- 
cesses [3, 2, 5] and the partition function in directed polymer models [2, 1, 11]. All 
these models belong to the so-called KPZ universality class, which is associated with 
unusual fluctuations of order t^^^ at time i on a spatial scale of i^/^. We refer the 
reader to the reviews [4, 12] for more details and background on the KPZ equation and 
universality class. 

As stated the KPZ equation (1.1) is ill-posed due to the non-linear term. To make 
sense of it we follow the approach of Bertini and Giacomin [3]. Observe that if we 
let Z{t,x) = exp{—T-L{t,x)) then, formally, Z{t,x) solves the (linear) stochastic heat 
equation with multiplicative noise 

(1.2) dtZ = ^d^Z - ZW. 

Therefore we simply deflne the solutions T-L{t,x) of (1.1) via the Hopf-Cole transforma- 
tion 



(1.3) 



^(t,x) = -\og{Z{t,x)), 
1 
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where Z{t, x) is the (weh-defined) solution of the stochastic PDE (1.2). In a remarkable 
recent development, M. Hairer [8] has proposed a way to make sense of the KPZ equation 
directly. The resulting solutions coincide with the Hopf-Cole solutions. 

One of the most interesting properties of the KPZ equation is the preservation of 
Brownian initial data. In particular, if one starts the equation with a standard two- 
sided Brownian motion, one sees at time t a new Brownian motion with the same 
diffusivity, but with a (random) height shift (the new Brownian motion will of course 
be coupled to the starting one in a highly non-trivial way). Furthermore, any initial data, 
however smooth, will immediately become locally Brownian. This can be understood 
in many ways. One is that one expects the local quadratic variation to be the same as 
that of the equilibrium solutions, for any positive time, for arbitrarily nice initial data. 
Another is that one expects that the solution at time t can be written as a standard 
two-sided Brownian motion B{x) plus a more regular object. One would naturally like 
to take this Brownian motion B{x) to be the solution of the equation starting from a 
two-sided Brownian motion. In other words, one would like to couple all solutions to 
the equilibrium one. For a large class of initial data, Hairer [8] has shown that the 

3 

solution can be written as a Brownian motion plus a function in Ca". Unfortunately, 
the Brownian motion used is a solution of the Langevin equation obtained by linearizing 
KPZ, as opposed to the equilibrium solution of KPZ itself. Moreover, not all initial data 
can be handled by the methods of [8] because they require that certain auxiliary objects 
be integrable against heat kernels in space and time. The singularity at time rules 
out one of the most important cases, which is the narrow wedge initial data. 

Our main interest will be this last case: the initial condition for KPZ given by starting 
the stochastic heat equation (1.2) with initial data 

(1.4) Z{0,x) = 6,=o. 

This defines T-L{t, x) for every t > via (1.3), but one should not think in terms of the 
initial data for T-L, since the delta function does not have a well-defined logarithm. This 
narrow wedge initial data is very basic. For example, it is the one that approximates 
the free energy of point-to-point polymers. 

For this initial data, if we define At{x) by 

n{t,x) = |! +log(\/2^ + ^ _ 2-1/3^1/3^^(2-1/3^-2/3^) 

then, properly rescaled, At{x) converges to a Gaussian process as t — )■ and it is 
conjectured that it converges to the Airy2 process as t — )■ oo (see Conjecture 1.5 in [2]). 
At{x) is therefore referred to as the crossover Airy2 process, and interpolates between 
the KPZ and Edwards- Wilkinson [6] universality classes, the last one associated with 
the stochastic heat equation with additive noise and hence Gaussian statistics. 

We will denote by T-L°'^{t, x) the solution of the KPZ equation (1.1) started with initial 
condition 

(1.5) ■H'''^iO,x)=B{x), 

where B{x) is a two-sided standard Brownian motion. We recall that this initial condi- 
tion is such that, for each fixed t > 0, T-l^^{t, x) — T-l°'^{t, 0) is itself a two-sided standard 
Brownian motion (in space), see Proposition B.2 in [3]. In this equation we use the same 
white noise as in the earlier solution starting with Dirac mass (1.4). The two solutions 
exist, are unique, and are coupled for all time. 

Our main result is the following: 
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Theorem 1. Fix t > and let T-L{t,x) and 'H''"^{t,x) be the Hopf-Cole solutions of the 
KPZ equation (1.1) with respect to the same white noise and with initial conditions given 
by (1.4) and (1.5). Then 'H{t,x) — (T-L^'^{t,x) — W^{t,Q)^ is a finite variation process. 

Our initial data is special, but it is in some sense the furthest possible from equilib- 
rium, and it has the benefit of a surprisingly simple proof. 

We remark that if 'H(t, x) is started with initial condition given by a two-sided Brow- 
nian motion plus a Lipschitz function then it is easy to show using our coupling method 
that it remains a two-sided Brownian motion plus a Lipschitz function for all t > 0. 
This follows from the results of Hairer [8], but the proof there is much more involved. 
The precise statement is given next, its short proof uses the same ideas as the proof of 
Theorem 1 and is given in Section 2. 

Theorem 2. Let 'H'^'^{t, x) andT-L{t,x) be the Hopf-Cole solutions of the KPZ equation 
(1.1) with respect to the same white noise and with initial conditions given respectively 
by (1.5) and 

n{0,x) = B{x) + ip{x), 

where B{x) is the same two-sided standard Brownian motion as in (1.5) and (p is a 
Lipschitz function. Then, for every fixed t> 0, n{t,x)-{W'i{t,x)-W'^{t,0)) is almost 
surely a Lipschitz function (with the same Lipschitz constant as <p). In particular, the 
law of7i{t,x) in a finite interval has finite relative entropy with respect to the law of 
B{x) in that interval. 



2. Proofs 



The proofs of Theorems 1 and 2 rely on considering the weakly asymmetric simple 
exclusion process, which provides a microscopic model for the KPZ process. The simple 
exclusion process with parameters p,q ^ [0, 1] (such that p + g = 1) is a {0, l}^-valued 
continuous time Markov process, where I's are thought of as particles and O's as holes. 
The dynamics of the process are as follows: each particle has an independent exponential 
clock with parameter 1; when the clock rings, the particle attempts a jump, trying to 
go one step to the right with probability p and one step to the left with probability q; 
if there is a particle at the chosen destination, the jump is supressed and the clock is 
reset. We refer the reader to [10] for a rigorous construction of this process. We will be 
interested in the case q > p, known as the asymmetric simple exclusion process (ASEP). 
More precisely, we will be interested in the weakly asymmetric simple exclusion process 
(WASEP), where we introduce a parameter in the model and let the asymmetry q — p 
go to with the parameter. 

Given any configuration rj G {0, 1}^ for the exclusion process we will denote by 
rj G { — 1, 1} the configuration given by rf^x) = 2r]{x) — 1 for each x G Z. To any simple 
exclusion process r]t we can associate the height function h{t, •) : M — > % in the following 
manner: 



f2iV(t)+ ifa:>0, 

Q<y<x 



(2.1) /i(t,x) 



2N(t) if X = 0, 

2iV(t)- ^ ^t{y) ifx<0, 

x<y<X) 

where N{t') is the net number of particles which crossed from the site 1 to the site 
up to time t. It is straightforward to check that the simple exclusion process can be 
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recovered from the height function by 

(2.2) ^t(x) = h{t, x) - h{t, X - 1). 

Next we introduce the scahng parameter e > 0, which should be thought of as going 
to 0, and consider WASEP with asymmetry e-*^/^, that is, 

We will denote by ryf the resulting WASEP, which we start with the step initial condition 

(2.3) r/o(a^) = lx>o- 

To r/f we associate the height function h^(t,x) via (2.1). Our main tool will be the 
convergence of a suitably rescaled version of to the solution of the KPZ equation 
with initial condition (1.4). 

The convergence of the height function was proved by Amir, Corwin, and Quastel [2] 
by performing a microscopic Hopf-Cole transform analogous to (1.3), an idea introduced 
originally by Gartner [7] and further developed in [3]. Let 

Is = \e-''\ Xs = \ Iog(f ) = + + o{e'l\ 

Ve=p + q-2^= ieV2 + 1,3/2 ^ o(e5/2). 

The Hopf- Cole transformed height function is given by 

(2.5) Ze{t,x) = 7£exp(^ - Xehe{e~'^t,£~^x) + Vety 

Observe that, with this definition, Z^{0,x) — )■ 6x=o as e — )• as discussed in Section 1.2 
of [2]. 

We regard the process Ze{t,x) as taking values in the space Z)([0, oo), Du(M)), where 
-Du(M) refers to right-continuous paths with left limits with the topology of uniform con- 
vergence on compact sets, which we endow with the Skorohod topology. The following 
result corresponds to Theorem 1.14 of [2]: 

Theorem 3. The family of processes (Zg)^^^ converges in distribution as e ^ Q in 
D([0, oo), /^^(M)) to the C{[Q,<x),C{K)) -valued process Z given by the solution of the 
stochastic heat equation (1.2) with initial condition Z(0,x) = 6x=o- 

We recall that the exclusion process is attractive, which for our purposes means that 
two copies r]^ and rj^ of the process with initial conditions 7?g < r]Q (which just means 
fJoi^) — ^o(^) ^o'^ ^) be coupled in such a way that Tjj < rjl for all t > 0. We will 
refer to this coupling as the basic coupling and refer the reader to [10] for more details. 

We will denote by rff^ a copy of WASEP in equilibrium, started with a product 
measure with density ^, and by h'^'^ the associated height function. To prove Theorem 
1 we will couple r/f with rfl^ using the basic coupling. The key result will be an estimate 
on the number of discrepancies between the two processes at time £~'^t in a window of 
size 0(e~^), Proposition 2.1 below. 

Let ry™"^ and ry™^^ denote copies of WASEP started with initial conditions 

%"^''^ = ^?oAC and r^r'' = vh{x)yrf^, 

where the minimum and maximum are meant sitewise. Observe that ry™™ corresponds 
to starting with no particles on the negative half-line and a product measure of density 
^ on the positive half-line, while 7]^'"^^ corresponds to starting with a product measure 
of density ^ on the negative half-line and all sites occupied on the positive half-line. We 
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will denote by /i™™ and /i™*^^ the height functions associated respectively to these two 
processes. 

Let Z^^^ and be the Hole-Copf transformed height functions associated to 

the corresponding initial conditions, which are defined in the same way as in (2.5) 
with the scaling (2.4) except that 7^ = 1. The proof in [2] of Theorem 3 can be adapted 
without difficulty (see [5] for the details) to show that that Z™'" and Z™^^ converge in 
distribution in D{[0, 00), Du{M.)) respectively to the solutions Z™"(t,x) and Z™'*'^(t,x) 
of the stochastic heat equation (1.2) with initial data 2™™(0, x) = exp(— B{x)^lx>o 
and -H™^^(0, x) = exp(— B{—x)^lx<:o, where B{x) is a standard one-sided Brownian 
motion. We define W^'''{t,x) = -\og{Z'^^''{t,x)) and W^''''{t,x) = - log(Z™^^(t, x)). 

Given any of the height functions h with the different initial conditions we are con- 
sidering, we will denote by its rescaled version 

he{t,x) = e^/^/i(e"^t,e"^x). 

Proposition 2.1. Assume r/| is started with the step initial condition (2.3) and fix 

a < b and t > 0. Then, under the basic coupling, 

a;e[a£;-i,6£-i]nZ 

< \ [/.--(t, b) - h^^it, a)] - \ \hf-{t, b) - hf%t, a)] 

almost surely. 

Proof. We construct the four processes Tyf , r]^^, r^™™ and r^™*™ together under the basic 
coupling, so attractiveness implies that 

for all t > 0. Using this we get by (2.2) that 

\ E \{hr^e-H,x) - hr-ie-\x- 1)) 

hf'^ie-H, x) - hf'^ie-H, x-1) 



2 



= ^ [^r^(*, b) - hr^it, a)] - [hf-it, b) - hf-it, a)] . 

Multiplying by e^/^ we obtain the desired bound. □ 

Proof of Theorem 1. Fix a finite interval / = [a,b] and let TV/(/) denote the total 
variation of / in I: 

n 

TVK/)= sup J^|/(x,)-/(x,_i)|. 
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Then clearly 

SO by Proposition 2.1 we get 

TYj{h,{t, .) - h'^.^t, •)) < [hr-'it, b) - h^^it, a)] - [JiTHt, b) - hf%t, a)] . 
On the other hand 

he{t,x) - hl\t,x) = -e''/^\-^[\og{ZS,x)) - log(Z|'i(t,x))] +e-V2A-ilog(7,). 

Thus by Theorem 3 and Theorem 2.3 in [3] we get that x) — }ip^{t, x)—e~'^/'^X~^ log(7£^ 
converges in distribution to T-L{t, x) — Ti'^'^^t, x) on the interval /. Note that this requires 
a very minor extension of the results of [2] and [3], namely that the processes /ig and 
hf^ built from exclusion processes running with the same background Poisson processes 
converge jointly to H and H^^. There are no issues involved in extending Theorem 3 
and Theorem 2.3 in [3] to this situation and therefore we omit the details. 

By the lower semicontinuity of TV/ we deduce that 

¥{TYj{n{t,-)- [H''^t,-)-n''''{t,o)]) >K) =¥{i:Yi{n(t,-)-'W\t,-)) > k) 

< limsupFfTV,f/i,(t,-) - }f^{t, •) - £-1/2^-1 iog(7j) > k] 



= limsupP(^TV/(^/ie(t,-) - }f^{t,-)^ > K 

< limsupP((/i^^"(t,6) - /ir"(^,a)) " (hf^it^b) - hf''{t,a)) > K 

To finish the proof of Theorem 1 we observe that the quantity inside the last probability 
above equals 

-e^/^A^i [log(Z--(t, b)) - log(Zf-(t, a)) - \og{Zr\t, b)) + log(Z--(t, a))] . 

Note the key point that the additive constants in (2.5) cancel. Using the convergence of 
and Z'^^^ to solutions of the stochastic heat equation discussed before Proposition 
2.1, the above converges in distribution to 6) - 'H"'"(t,a)] - 6) - 

a)] . Since the last random variable is finite we deduce that 

^lini P(TV/(?^(t,-) - -^°'^(t,0)]) > K) =0. □ 

Proof of Theorem 2. Let r]t be a copy of WASEP started with product measure with 
density profile given by 

P(770(X + 1) = 1) = i + ie-V2 (^(^^) _ ^(^(^ _ 1))) . 

Then Theorem 2.3 of [3] implies that hs{t,x) converges in distribution as e — t- in 
£'([0, oo), £)„(M)) to T-L^t, x) (with initial condition as in the statement of the theorem). 
Now denote by M the Lipschitz constant of ip and let hf and denote the rescaled 
height functions corresponding to WASEP started respectively with product measures 
of densities |(1 + e^/'^M) and ^(1 — e^/^M). Coupling the initial conditions in the 
natural way and using the basic coupling and attractiveness, it is clear that h~{t,x) < 
he{t, x) < hf{t, x) for all t > 0. On the other hand, since product measures are invariant 
for WASEP, hf{t , x) converges in distribution to T-L'^'^{t, x) it Mx, and as before this 
convergence can be achieved jointly for h^, hf and h~ . Therefore, given any a < b, 
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H^^itjx) — Mx < 'H{t,x) < VP'^{t,x) + Mx almost surely for every x e [a, 6], and the 
result follows. □ 
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